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ABSTRACT 

Aims. We aim at studying the cluster environment surrounding the UCH n region in IRAS 20293+3952, a region in the first stages of 
formation of a cluster around a high-mass star. 

Methods. BIMA and VLA were used to observe the 3 mm continuum, N 2 H + (1-0), NH 3 (1,1), NH 3 (2,2), and CH 3 OH (2-1) emission 
of the surroundings of the UCH n region. We studied the kinematics of the region and computed the rotational temperature and column 
density maps by fitting the hyperfine structure of N 2 H + and NH 3 . 

Results. The dense gas traced by N 2 H + and NH 3 shows two different clouds, a main cloud to the east of the UCH n region, of ~0.5 pc 
and ~250 Mq, and a western cloud, of ~0.15 pc and ~30 Mq. The dust emission reveals two strong components in the northern 
side of the main cloud, BIMA 1 and BIMA 2, associated with YSOs driving molecular outflows, and two fainter components in the 
southern side, BIMA 3 and BIMA 4, with no signs of star forming activity. Regarding the CH 3 OH, we found strong emission in a 
fork-like structure associated with outflow B, as well as emission associated with outflow A. The YSOs associated with the dense gas 
seem to have a diversity of age and properties. The rotational temperature is higher in the northern side of the main cloud, around 
22 K, where there are most of the YSOs, than in the southern side, around 16 K. There is strong chemical differentiation in the region, 
since we determined low values of the NH 3 /N 2 H + ratio, ~50, associated with YSOs in the north of the main cloud, and high values, 
up to 300, associated with cores with no detected YSOs, in the south of the main cloud. Such a chemical differentiation is likely due 
to abundance/depletion effects. Finally, interaction between the different sources in the region is important. First, the UCHii region 
is interacting with the main cloud, heating it and enhancing the CN(l-O) emission. Second, outflow A seems to be excavating a 
cavity and heating its walls. Third, outflow B is interacting with the BIMA 4 core, likely producing the deflection of the outflow and 
illuminating a clump located ~0.2 pc to the northeast of the shock. 

Conclusions. The star formation process in IRAS 20293+3952 is not obviously associated with interactions, but seems to take place 
where density is highest. 
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1. Introduction 

It is generally accepted that the formation of massive stars takes 
place not isolated but simultaneously with the formation of a 
cluster. Clusters of infrared sources have been studied at op- 
tical and near-infrared wavelengths around intermediate/high- 
mass stars (e.g., McCaughrean & Stauffer 1994; Hillenbrand 
1995; Hillenbrand & Hartmann 1998; Testi et al. 1998, 1999, 
2000). One could surmise on the structure and the evolution- 
ary state of the molecular cloud when such clusters were being 
born around an intermediate/high-mass protostar. In the standard 
theory of star formation, the central massive object, deeply em- 
bedded in large amounts of gas and dust, starts to radiate at UV 
wavelengths at very early stages, before finishing the accretion 
of all its final mass (e. g., Bernasconi & Maeder 1996). The UV 
photons ionize the surrounding medium, and the high-mass star 
developes an UCHii region around it (see, eg. Garay & Lizano 
1999), which expands and pushes forward the gas and dust sur- 
rounding it until the parental cloud is completely disrupted and a 
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cluster of infrared/optical sources emerges (e. g., Hester & Desch 
2005). The study of the clumpy medium around massive proto- 
stars helps not only in the understanding of the formation of the 
massive star itself, but also provides a direct way to study the 
clustered mode of star formation, in which most stars are thought 
to form. 

In order to properly characterize the medium around a mas- 
sive protostar, it is necessary to be sensitive to small (~2000 AU) 
and low-mass (55 1 Mq) clumps of gas and dust. For this, in- 
terferometric observations toward relatively nearby (< 3 kpc) 
massive star-forming regions are required, and are now possible 
with sufficient sensitivity. Here we report on observations of one 
intermediate/high-mass star-forming region, IRAS 20293+3952, 
located at 2.0 kpc of distance (Beuther et al. 2004b) in the south- 
western edge of the Cygnus OB2 association, and with a lumi- 
nosity of 6300 Lq. This region has been included in the sam- 
ple of high-mass protostar candidates of Sridharan et al. (2002). 
The IRAS source is associated with the only centimeter source 
detected in the region, which is likely tracing an UCH n region 
(Beuther et al. 2002b). A study of the H2 emission in the re- 
gion shows two near-infrared stars (IRS 1 and IRS 2) associated 
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with the IRAS source, and a circular ring of H2 emission sur- 
rounding IRS 1 (Kumar et al. 2002). Beuther et al. (2002a) ob- 
served the region with the IRAM 30 m telescope and find some 
substructure at 1 .2 mm around the UCH n region. The strongest 
millimeter peak, very close to the position of a H2O maser, is 
located ~15" north-east of the IRAS source, and observations 
with high angular resolution by Beuther et al. (2004a) reveal 
a compact and strong millimeter source, mml, associated with 
the H2O maser. Two other fainter compact millimeter sources, 
mm2 and mm3, are located 10" to the east of the UCH n region. 
Beuther et al. (2004a), from CO (2-1) and SiO(2-l) observa- 
tions, suggest the presence of four molecular outflows in the re- 
gion, with two of them, outflows A and B, associated with two 
chains of H2 knots (Kumar et al. 2002). Subsequent observations 
with the PdBI were carried out at 2.6 and 1.3 mm by Beuther 
et al. (2004b), who find CN emission close to mml, mm2 and 
mm3. The different millimeter sources detected and the pres- 
ence of multiple outflows around the UCH n region makes this 
region a good choice to study star formation and interaction in a 
clustered environment. 

In this paper we report on BIMA and VLA observations of 
the continuum at 3 mm and the dense gas traced by N2H + (1-0), 
NH 3 (1,1), and NH 3 (2,2) toward IRAS 20293+3952, together 
with observations of several CH30H(2-1) transitions. In §2 
we summarize the observations and the reduction process, in 
§ 3 we show the main results for the continuum and molecu- 
lar line emission, in § 4 we analyze the line emission and show 
the method used to derive the rotational temperature and col- 
umn density maps. Finally, in § 5 we discuss the results obtained, 
mainly the properties of the dense gas surrounding the UCH n 
region, the different sources identified in the region, and the in- 
teraction between them. 

2. Observations 

2.1. BIMA 

The N 2 H + (l-0) and CH 3 OH(2-l) lines and continuum at 
95 GHz were observed towards the IRAS 20293+3952 region 
with the BIMA arra)Q at Hat Creek. The observations were car- 
ried out on 2003 September 28 in the C configuration, and on 
2004 March 24, in the B configuration, with 10 antennas in use. 
The phase center of the observations was set at a(72000) = 
20 h 31 m 12 s .70; 5(72000) = +40°03'13'/4 (position of the mil- 
limeter peak detected with the IRAM 30 m telescope by Beuther 
et al. 2002a). The full range of projected baselines, including 
both configurations, was 9.4-240 m. The FWHM of the primary 
beam at the frequency of the observations was ~120". Typical 
system temperatures were ~500 K for Sep 28 and 500-1500 K 
for Mar 24. 

The digital correlator was configured to observe simultane- 
ously the continuum emission, the N2H + (1-0) group of hyper- 
fine transitions (93.17631 GHz, in the lower sideband), and the 
CH3OH 2 1 — 1 1 E, 2()-lo E, 2o-lo A + , and 2 1 — 1 1 E transitions 
(96.74142 GHz, in the upper sideband). For the continuum, we 
used a bandwidth of 600 MHz in each sideband, and for the lines 
we used a bandwidth of 25 MHz with 256 channels of 100 kHz 
width, providing a spectral resolution of 0.31 km s _1 . 

Phase calibration was performed with QSO 2013+370, with 
typical rms in the phases of 12° and 43° for C and B config- 
urations, respectively. The absolute position accuracy was esti- 
mated to be around a/5. We also used QSO 2013+370 for flux 
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calibration, and the error in the flux density scale was assumed to 
be ~30%. Data were calibrated and imaged using standard pro- 
cedures in MIRIAD (Sault et al. 1995). We combined the data 
from B and C configurations. In order to improve the angular 
resolution of the continuum emission, we weighted the visibility 
data with a robust parameter of +1. For the line emission, which 
is more extended than the continuum, we used natural weight- 
ing. The resulting synthesized beams and the final rms noises 
are listed in TableQ] 

2.2. VLA 

Observations of the (7, K) = (1,1) and (2,2) inversion transi- 
tions of the ammonia molecule were carried out with the Very 
Large Array (VLA) of the NRAC0 in the D configuration on 
2000 September 3. The phase center was set to a(72000) = 
20 h 31 m 10 s .70; 5(72000) = +40°03'1077 (catalogue position of 
the IRAS source). The FWHM of the primary beam at the fre- 
quency of observation was ~1 10", and the range of projected 
baselines was 37.9-708 m. The absolute flux calibration was 
performed by using 3C48, adopting a flux density at 1.3 cm of 
1.05 Jy. The phase calibrator was QSO 2013+370, with a 1.3 cm 
bootstrapped flux density of 3.84 Jy, and 3C 84 was used as the 
bandpass calibrator. 

The NH 3 (1,1) and NH3 (2,2) lines were observed simulta- 
neously in the 4 IF correlator mode of the VLA (with 2 po- 
larizations for each line), providing 63 channels with a spec- 
tral resolution of 0.62 km s across a bandwidth of 3.13 MHz, 
plus an additional continuum channel containing the central 75% 
of the total bandwidth. The bandwidth was centered at the sys- 
temic velocity vlsr = 6.3 km s (Sridharan et al. 2002) for the 
NH 3 (1,1) line, and at v LSR = 1 1.3 km s" 1 for the NH 3 (2,2) line 
(to cover the main and one of the satellite components). Data 
were calibrated and imaged using standard procedures of AIPS. 
The cleaned channel maps were made using natural weighting 
of the visibility data. Table Q] summarizes the parameters of the 
observations. 



3. Results 

3.1. Continuum emission 

Figure [T] shows the continuum map at 3.15 mm overlaid on 
the continuum infrared emission at 2.12 /urn from Kumar et al. 
(2002). The millimeter continuum emission has two strong com- 
ponents, BIMA 1 and BIMA 2, both spatially resolved. While 
BIMA 1 is elongated in a northeast-southwest direction (roughly 
PA.^45°), BIMA 2 is somewhat flattened in the east-west di- 
rection, and is surrounded by an elongated structure of ~30" 
(0.3 pc), which is tracing a dust ridge. There is an extended 
clump 10" to the southeast of the dust ridge, which we label 
BIMA 3. Note that in this southeastern region there are very few 
infrared sources. At 10" to the east of the ridge, there is an- 
other faint feature, BIMA 4, extending for 10" in the northwest- 
southeast direction. Continuum emission at a level of 5<x is de- 
tected toward the position of the UCH 11 region, and toward 
IRS 2. Beuther et al. (2004b) observed the same region at 2.6 mm 
with the PdBI at higher angular resolution ( 1'.'5 x l'/2), but with a 
smaller primary beam (44"). The authors find a very strong com- 
pact source, mml, toward the position of BIMA 1 and elongated 
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PA. 


Spec, resol. 


Rms a 


Observation 


Telescope 


Config. 


(arcsec) 


O 


(km s" 1 ) 


(mJy beam - ') 


Continuum 


BIMA 


B+C 


5.8x5.6 


-6.4 




0.7 


CH 3 OH (2-1) 


BIMA 


B+C 


6.1 x 5.7 


-6.2 


0.31 


80 


N 2 H + (l-0) 


BIMA 


B+C 


6.2x5.9 


-4.1 


0.31 


70 


NH 3 (U) 


VLA 


D 


6.9x3.0 


71.5 


0.62 


3 


NH 3 (2,2) 


VLA 


D 


6.6x3.1 


71.4 


0.62 


3 



a Rms noise per channel in the case of line emission. 




Sight Ascension (J2000) 

Fig. 1. Contours: 3.15 mm continuum emission from IRAS 20293+3952 obtained with BIMA in the B and C configurations using 
an intermediate weight between uniform and natural (robust = +l).The contour levels are -3, 3, 5, 7, 9, 12, 15, 18, 21, and 24 times 
the rms noise of the map, 0.67 mJy beam -1 . The synthesized beam, shown in the bottom right corner, is 5'.'8 x 5'.'6, at P.A. = -6°4. 
White contours correspond to the centimeter emission at 3.6 cm, and trace the UCHn region (Beuther et al. 2004a). Contours are 
-3, 3, 6, and 9 times 70 fiiy beam -1 . Grey scale: Continuum image at 2.12 /mi from Kumar et al. (2002). Main clumps of dust 
emission are labeled as BIMA 1 to BIMA 4. White dots indicate the compact millimeter sources detected by Beuther et al. (2004b) 
at 2.6 mm with the PdBI. These are, from west to east, mm2, mm3, and mml. White crosses mark the sources from 2MASS with 
infrared excess, and black tilted crosses mark the position of the Northern Warm Spot (NWS) and the Southern Warm Spot (SWS). 



in the same northeast-southwest direction as BIMA 1 . Besides 
this strong component, a weak subcomponent was resolved in 
the direction of the elongation (to the southwest of mml), and 
was labeled mm la. Beuther et al. (2004b) also find two compact 
millimeter sources associated with BIMA 2, mm2 and mm3 (see 
Fig.Q]for the positions of these compact millimeter sources). 

In Table [2] we show the position, peak intensity, flux density 
and masses associated with BIMA 1, BIMA 2, the entire dust 
ridge (including BIMA 2), BIMA 3, BIMA 4, and IRS 2. We 
discuss these results in section |4~T1 



3.2. CHiOH 

Figure [2^ shows the zero-order moment map of the CH3OH 
emission, integrated from -1.5 to 21 km s -1 , including the two 
strongest transitions 2o-lo A + , and 2_i-l_i E. Figure [3] shows 
the CH3OH spectrum at selected positions. The channel maps 



corresponding to the lines 2o-lo A + and 2_i — 1 1 E are shown in 
Fig.E] 

The strongest emission of CH3OH is found to the southeast 
of BIMA 1, elongated in the southeast-northwest direction, and 
covering a spatial extension of ~55" (0.6 pc). The emission has 
a fork-like structure (this is well observed in the 7.7 km s -1 
velocity channel of Fig. HJ, and extends through several chan- 
nels. This structure is associated with high-velocity SiO(2-l) 
and CO (2-1) emission of outflow B (Beuther et al. 2004a). The 
position-velocity (p-v) plot along outflow B is shown in Fig. 
The most prominent feature in the plot is a blue wing, ~6 km s -1 
wide, spanning from 0" to 15" offset positions. At negative off- 
sets the CH3OH emission is dominated by a weak red wing 
(this is better seen for the 2 1 — 1 1 E line). Farther away than 
~25" (which corresponds to the southeasternmost clump of out- 
flow B), the emission has contributions from both redshifted and 
blueshifted emission. This change of behaviour spatially coin- 
cides with BIMA 4, located in the p-v plot at an offset position 
of 20". 
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Fig. 2. a) CH 3 OH zero-order moment for the 2 -lo A + , and 2 1— 1 1 E lines toward IRAS 20293+3952. Contours start at 1%, 
increasing in steps of 12% of the peak intensity, 10.7 Jy beam 1 km s _I . b) N2H + zero-order moment integrated for all the hyperfine 
components of the (1-0) transition. Contours start at 2%, increasing in steps of 12% of the peak intensity, 21.9 Jy beam -1 km s _1 . 

c) NH3 (1,1) zero-order moment. Contours start at 5%, increasing in steps of 15% of the peak intensity, 0.240 Jy beam -1 km s _1 . 

d) NH3 (2,2) zero-order moment. Contours start at 7%, increasing in steps of 15% of the peak intensity, 0.0892 Jy beam -1 km s _1 . 
In all panels, the first level is about 3 times the rms noise of the map. Symbols are the same as in Fig. Q] with the star marking 
the peak of the centimeter emission, and the triangles marking the positions of BIMA 3 and BIMA 4 cores. MEC and MWC stand 
for methanol eastern clump and methanol western clump, respectively. The synthesized beams for each transition are shown in the 
bottom right corner, and are listed in TableQ] The straight lines mark the direction of outflows A and B (Beuther et al. 2004a), with 
the light grey line corresponding to the redshifted lobe. The dashed curve indicates the primary beam of BIMA (panels a and b) and 
VLA (panels c and d) observations. 
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Table 2. Parameters of the continuum sources at 3.15 mm 
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Peak Position 




Flux" 


Assumed 7d b 


Mass 




o-(J2000) 


c5(J2000) 


(mJy beam 1 ) 


(mJy) 


(K) 


(Mq) 


BIMA1 


20:31:12.767 


40:03:22.65 


17.4 


28 


34 


6.3 


BIMA2 


20:31:12.062 


40:03:12.22 


12.1 


23 


34 


5.2 


BIMA 2+ridge 








42 


25 


14 


BIMA 3 


20:31:13.769 


40:03:03.12 


5.4 


7.9 


17 


3.9 


BIMA 4 


20:31:14.558 


40:03:06.95 


4.2 


7.9 


17 


3.9 


IRS 2 


20:31:10.303 


40:03:17.10 


3.5 


3.9 


50 


0.60 



a Flux density inside the 5<r contour level for BIMA 1, BIMA 2, and BIMA 2+ridge, and the 3cr level for BIMA 3, BIMA 4, and IRS 2. 

b r d is estimated by correcting the rotational temperature derived from NH 3 (see § 14.21 ) to kinetic temperature, following the expression of 

Tafalla et al. (2004). For IRS 2 we assumed T A ~ 50 K. 
c Masses of gas and dust derived assuming a dust emissivity index /3 = 1 (Beuther et al. 2004b). The uncertainty in the masses due to the opacity 

law and the dust emissivity index is estimated to be a factor of four. 



A clump to the east of outflow B that can be seen in the inte- 
grated intensity map shows narrow lines and appears only at sys- 
temic velocities (see Figs. [3] and @]). We label this clump, which 
has not been previously detected, 'methanol eastern clump'. 

In addition, there are two redshifted clumps, one associated 
with BIMA 1, and the other close to the UCH n region (see, e. g., 
velocity channels from 13 to 16 km s~'), which probably are 
part of the redshifted CO (2-1) lobe of outflow A (Beuther et al. 
2004a). We label the clump close to the UCH n region 'methanol 
western clump'. In Fig. [6] showing the p-v plot along outflow A, 
we find that there is a velocity gradient, from 4 km s _1 at 7", to 
8 km s _1 at -10", clearly seen in the two lines. This velocity gra- 
dient is centered around the zero offset position (corresponding 

to mml). The methanol western clump, seen at 35", is highly 

redshifted (line 2o-lo A + is at ~9 km s~') and shows broad red 
wings, spreading ~5 km s _1 . 

3.3. N 2 tT 

The zero-order moment map integrated for all the hyperfine tran- 
sitions is presented in Fig. [2b. Figure [3] shows the N2FU (1-0) 
spectra at selected positions. 

The integrated N2FU emission consists of a main cloud and 
a smaller and weaker cloud to the west of the main cloud (west- 
ern cloud). The size of the main cloud and the western cloud are 
-50" (0.5 pc), and -15" (0.15 pc), respectively. The four BIMA 
sources are located in the main cloud with BIMA 1 close to the 
N2H + emission peak. N2H + emission is marginally detected to- 
wards the UCH 11 region. 

The hyperfine F\ = 0-1 was used for the analysis of the 
kinematics of the N2H + emission because it is not blended with 
the other hyperfines. In Figs.|7]and[8]we show the channel maps 
and the first and second-order moments for this hyperfine. The 
channel with maximum intensity was found at the systemic ve- 
locity (6.3 km s _1 ). 

The first-order moment map shows a small velocity gradi- 
ent in the main cloud with increasing velocities roughly from 
the southwest to the northeast. In addition, there is a filament 
spatially coincident with outflow B (clearly visible in the 7.1 to 
7.7 km s _I channel maps of Fig- EJ, which has very broad lines 
associated, as seen in the second-order moment map (Fig. [8}?). 
The p-v plot made across outflow B and within the main cloud 
(Fig-Ep) shows the aforementioned velocity gradient in the main 
cloud and the line broadening toward outflow B (at offset posi- 
tion zero). The p-v plot along outflow B (Fig. |5};) shows clear 



line broadening all along the outflow, but without the wing emis- 
sion component found in the CH3OH emission. 

The western cloud appears redshifted, at > 7 km s _1 , and has 
narrow lines (Fig.[7]l. At these velocities, there is weak emission 
connecting the western cloud with outflow B. Figure [9] shows 
the p-v plot along a cut at P.A.=110°, picking up partially the 
outflow and the western cloud. From this plot, it seems that the 
western cloud extends ~V (0.6 pc) toward the east, intersecting 
the main cloud. 

3.4. NH 3 

The zero-order moment map of the NH3 (1,1) emission (inte- 
grated including the main line and the inner satellites) is shown 
in Fig. |2};. The overall structure of NH3 (1,1) resembles roughly 
the emission of the dust and the emission of N2H + . NH3 (1,1) 
emission is also associated with the main cloud and the west- 
ern cloud. There is another NH3 (1,1) clump to the northwest 
(not shown in Fig. |2), which is very weak in N2H + (in part due 
to the primary beam attenuation). A difference between N2H + 
and NH3 (1,1) emission is that whereas N2H + peaks close to the 
millimeter continuum sources, the strongest NH3 (1,1) emission 
in the main cloud lies in the southeastern part. In addition, the 
strongest NH3 (1,1) emission over the entire region comes from 
the western cloud, not from the main cloud. NH3 (1,1) emission 
toward the UCH 11 region is weak. 

Figure [2ji shows the zero-order moment map of the 
NH3 (2,2) emission. NH3 (2,2) emission resembles closely that 
of the NH3 (1,1). However, the main difference with NH3 (1,1) 
is that the strongest emission of the NH3 (2,2) in the entire re- 
gion is found very close to BIMA 1 . The NH3 (2,2) emission 
in the western edge of the main cloud shows an extension that 
encompasses the position of the UCH 11 region. 

Note that the dense gas traced by NH3 and N2H + is morpho- 
logically very different from the gas traced by CH3OH. 

4. Analysis 

4.1. Dust 

We compared the position of BIMA 1 with the millimeter source 
found by Beuther et al. (2004b) at 3.5 mm at similar angular 
resolution. The peak of BIMA 1 is shifted ~l" to the west of 
mml. The values of the peak intensities (Table [2]) are in very 
good agreement with Beuther et al. (2004a), and the flux den- 
sities from Table [2] are about 50% larger than those reported in 
Beuther et al. (2004a). The offset in positions and the larger flux 
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Fig. 3. Spectra toward seven positions of the IRAS 20293+3952 region for the four transitions studied in this paper, from left to 
right, NH3 (1,1), NH3 (2,2), N2H + (1-0), and CH3OH (2-1). The seven positions are labeled on the right panel of each row, and are, 
from top to bottom, BIMA 1, BIMA 3, BIMA 4, NWS (Northern Warm Spot), SWS (Southern Warm Spot), MEC (methanol eastern 
clump), and MWC (methanol western clump). The vertical scale is indicated for each transition in the bottom row. For CH3OH we 
show, in order of increasing velocity, the 2o-lo E, 2o-lo A + , and 2 1 — 1 1 E lines (velocities are referred to the 2q-1q A + line). For 
N2H + , velocities are referred to the F\ = 0-1 hyperfine. 



densities detected in this work could be produced by the differ- 
ent spatial filtering of the BIMA and PdBI arrays, with the BIMA 
array more sensitive to large-scale structures. 

The continuum emission at 3. 15 mm is probably due to ther- 
mal dust emission for most positions of the region, since free- 
free emission from ionized gas, traced by the continuum emis- 
sion at 3.6 cm (Fig.Q]), is detected only at the western edge of 
the main cloud, at the position of the UCH 11 region. However, 
Beuther et al. (2004b) find a spectral index between 1.3 and 
2.6 mm toward mmla (the subcomponent to the west of mml) of 
0.8, and suggest that the collimated ionized wind from outflow 
A could account for such a spectral index. The contribution of 



the possible free-free emission from mmla to BIMA 1 is small, 
given that the intensity of mml is five times larger than the in- 
tensity of mmla, and thus essentially all continuum emission at 
3 mm is due to thermal dust emission. 

To derive masses from the flux densities at 3. 15 mm, we cor- 
rected the rotational temperature, between 15 and 25 K (calcu- 
lated from NH 3 (1,1) and NH 3 (2,2), see §03]!, to kinetic tem- 
perature (Walmsley & Ungerechts 1983), by following Tafalla et 
al. (2004), yielding kinetic temperatures in the range 17-34 K. 
The masses (Table |2]i are 2 times larger than those obtained by 
Beuther et al. (2004b) from observations at 3.5 mm. The differ- 
ence arises from the fact that BIMA is detecting more flux than 
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Fig. 4. CH3OH channel maps toward IRAS 20293+3952 for the lines 2o-lo A + (around 6.3 km s~', the systemic velocity), and 
2_i-l_i E (around 13 km s _1 ), averaged every 4 channels. Contours are -16, -12, -8, -4, 4, 8, 12, 16, 20, 24, 28, and 32 times the 
rms of the map, 0.04 Jy beam 1 . Filled circles are the compact millimeter sources detected by Beuther et al. (2004b), and the star 
marks the position of the UCH n region. The synthesized beam is shown in the bottom right corner of each panel. 



PdBI observations, and from the different dust temperatures as- 
sumed. While Beuther et al. (2004a) assume a dust temperature 
(derived from a graybody fit to the spectral energy distribution 
in the entire region) of 56 K, we adopted the kinetic tempera- 
tures obtained from the NH3 observations, which allowed us to 
separate the contribution from each source, thanks to the high 
angular resolution. The lower temperatures adopted imply larger 
masses in order to produce the same flux density. 

4.2. Rotational temperature and column density maps 

We obtained N2H + and NH3 spectra for positions in a grid of 
2" x 2" in the main cloud and the western cloud, and fitted the 
hyperfine structure of each spectrum, or a single Gaussian for 
NH3 (2,2). We performed the fits only for those spectra with an 
intensity greater than 5cr for NH3 (1,1), and greater than 4cr for 
NH 3 (2,2). We set a higher cutoff for NH 3 (1,1) to ensure that 
not only the main line is detected but also the satellites. The 
results from the fits indicate that the entire main cloud is essen- 
tially optically thin for N2H + (ttot ^ 1.5), but optically thick for 
NH3 (1,1) (ttot ^ 12). In both cases, the highest opacities are 
reached in the southeastern side of the cloud, around BIMA 3. 

From the results of the fits to the NH 3 ( 1 , 1 ) and the NH 3 (2,2) 
spectra we computed the rotational temperature and NH3 col- 
umn density maps. We derived the rotational temperature fol- 
lowing the standard procedures (Ho & Townes 1983; Sepulveda 
1993), and assuming that the excitation temperature and the 



intrinsic line width were the same for both NH3(1,1) and 
NH3 (2,2). For the NH3 column density derivation, we followed 
the expression given in Table 1 of Anglada et al. (1995), and 
assumed that the filling factor was 1 . 

The map of the rotational temperature is presented in 
Fig.fTOh. The maximum value, 38 ± 15 K, is reached at the po- 
sition of the UCH 11 region. Interestingly, around the apparent 
dense gas cavity west to BIMA 1, NH3 shows high tempera- 
tures, with a maximum of 34 ± 9 K (the fits to the NH3 (1,1) and 
NH3 (2,2) at these positions are reasonably good, so this heating 
must be real). Toward the western cloud we find that temperature 
is essentially constant, around 16 K, slightly decreasing toward 
the center of the cloud. 

There are some local maxima of rotational temperature in 
the map. In order to test the significance of these local maxima, 
we fitted them with a 'gaussian + background level' model in 
a region a few times the synthesized beam size, with the gaus- 
sian width and position angle fixed and equal to the synthesized 
beam. We defined the temperature enhancement as the difference 
between the local maxima and the background level. The rms of 
the residual image in the fitted region is < 1.5 K, and the posi- 
tions of the gaussian peaks are within 072 of the local maxima 
positions. In Table[3]we list the positions of the significant local 
maxima, their temperature, and temperature enhancement. We 
considered that a local maxima is significant when the temper- 
ature enhancement is at least two times the uncertainty. There 
is a local maximum that is nearly coincident (within 0'.'5) with 
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Fig. 5. a) CH3OH position-velocity (p-v) plot along outflow B (P.A.=130°). Channel maps have been convolved with a beam of 
5"x2", with P. A. perpendicular to the direction of the cut. Contours start at 0.17 Jy beam -1 , and increase in steps of 0.17 Jy beam -1 . 
The bottom dashed line indicates the systemic velocity for line 2o-lo A + (taken as the reference line), at 6.3 km s _1 . The top dashed 
line indicates the systemic velocity for line 2 1 — 1 1 E, at 12.6 km s _1 . b) N2H + p-v plot for the F\ = 0-1 hyperfine across outflow 
B and the main cloud (PA. =40°). Positive positions are toward the northeast. Contours start at 0.2 Jy beam -1 , increasing in steps 
of 0.2 Jy beam 1 , c) N2H + p-v plot for the F\ = 0-1 hyperfine along outflow B (PA. =130°). The velocity scale is placed matching 
the velocity scale of the CH3OH emission. Contours start at 0.15 Jy beam 1 , increasing in steps of 0.15 Jy beam 1 . For N2H + , we 
did not convolve the channel maps. In all panels, the central position corresponds to the Northern Warm Spot (see Tabled, and is 
marked by a vertical dashed line in panel a. 



IRS 5, with a value of 24 ± 3 K. The warmest local maximum is 
located ~5" to the southeast of mml (labeled as 'Northern Warm 
Spot'), reaching 29 + 6 K, and is apparently associated with faint 
continuum emission at 2.12 /mi (see Fig. [T). There is another 
local maximum associated also with faint emission at 2.12 jjm, 
which is found 10" to the southwest of BIMA 3, with 20 + 2 
K (labeled as 'Southern Warm Spot'). Another local maximum 
is associated with the H2 knot c (Kumar et al. 2002), lying on 
the axis of outflow B. Finally, note that there is also some heat- 
ing about 20" to the southeast of the Northern Warm Spot, but 
we did not consider this heating in the table because it falls in 
the edge of the region where we fitted the spectra. Note that all 
the temperature enhancements in Table[3]have infrared emission 
associated, giving support to their significance. 

FigurefTUb shows the resulting column density map for NH3, 
corrected for the primary beam response. An obvious feature 
from the NH3 column density map is that the highest values are 
found to the southeast of the main cloud, where we found the 



lowest values in the rotational temperature map (Fig.fTOb). In the 
northern part of the cloud the column density is higher around 
the Northern Warm Spot. It is worth noting that the NH3 column 
density decreases toward IRS 5. Note also that in the western 
cloud the column density increases toward the center. 



We calculated the NJti + column density by following 
Benson et al. (1998), taking into account the opacity effects, and 
correcting for the primary beam attenuation. In the expression 
for the column density, the value for T ex was derived from the 
hyperfine fit made with CLASS, and assuming a filling factor of 
1. The resulting map is shown in Fig.fTOfc. Contrary to NH3, the 
map of the N2FL column density has the maximum value very 
close to the position of the Northern Warm Spot and BIMA 1, 
and not in the southern side of the main cloud. 
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Fig. 6. CH 3 OH p-v plot along outflow A, with P.A.=67°. Contours are -0.25, -0.15, 0.15, 0.25, 0.35, 0.45, and 0.52 Jy beam 1 . The 
velocity resolution is 0.30 km s _1 , and the channel maps have been convolved with a beam of 10" x 5", with a P.A. perpendicular 
to the direction of the cut, in order to recover the maximum emission in each position. The velocity range includes two lines: 2o-lo 
A + (at 6.3 km s _1 , bottom dashed line) and 2_i-l_i E (at 12.6 km s _1 , top dashed line). The zero position corresponds to mml, and 
positive values go to the northeast. 
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4.3. The NH 3 /N 2 H + column density ratio map 

In order to compare the NH3 emission with the N2H + emission, 
we convolved the NH3 and N 2 H + channel maps to obtain a final 
beam of 7" (the major axis of the NH3 and N 2 H + beams). We 
fitted the spectra in each position of a grid of 4" x 4" in the con- 
volved maps, and derived the column density for NH3 and N 2 H + 



following the same procedures outlined above, and correcting 
for the primary beam of each interferometer. We computed the 
NH3/N 2 H + column density ratio map, and the result is shown in 
Fig.fTOb. From the NH3/N 2 H + ratio map, one can see a clear gra- 
dient from the northwest of the main cloud, with a ratio around 



10 Palau et al.: Star Formation in IRAS 20293+3952 




RA (J2000) 

Fig. 8. a) Color scale: first-order moment map for the hyperfine F l - 0-1 line of N 2 H + (1-0) toward IRAS 20293+3952. b) Color 
scale: second-order moment map for the hyperfine F\ = 0-1 line of N2H + (1-0). In both figures black contours are the same as 
Fig.|2k, showing the CH30H(2-1) emission, with contours starting at 13%, and increasing in steps of 15% of the peak intensity. 
Color scales are in m s -1 . The synthesized beam is shown in the bottom right corner, and symbols are the same as in Fig. [2] Note 
that the second-order moment gives the velocity dispersion, and must be multiplied by the factor 2 V21n2 ^ 2.35 to obtain FWHM 
line widths. 



Table 3. Significant local temperature enhancements 



Position 


T m t 


Ar 1 


Possible 


a(J2000) (5(J2000) 


(K) 


(K) 


Counterparts 


20:31:13.26 40:02:58.9 


20 ±2 


7±2 


2.12 yum faint emission (SWS b ) 


20:31:13.23 40:03:19.5 


29 ±6 


12 ±6 


2.12 fim faint emission (NWS b ) 


20:31:13.44 40:03:13.4 


24 ±3 


7±3 


IRS 5 


20:31:14.96 40:03:03.4 


23 ±4 


9±4 


H2 knot c c 



a We assumed that each local maximum can be described by a 'gaussian + background level' model, with the gaussian width and position angle 
fixed and equal to the synthesized beam. The temperature enhancement, Ar, is defined as the difference between the local maximum and the 
background level. We considered as significant the temperature enhancements with Ar > 2cr, with a being the uncertainty. 

b NWS: Northern Warm Spot; SWS: Southern Warm Spot 

c from Kumar et al. (2002). 



50, to the southeast, where the ratio reaches values up to ~300. 
Such high values are also reached in the western cloud. 

5. Discussion 

5.1. General properties of the dense gas 
5.1.1. Rotational temperature 

The temperature distribution in the main cloud can be clearly 
separated into two parts: the northern side, with an average tem- 
perature of ~22 K, and the southern side, with an average tem- 
perature of ~16 K. It is interesting to note that almost all the 
YSOs in the region (§ !5.3b are associated with the northern side 
of the cloud, while in the cold southern side we found very few 
hints of active star formation. Thus, the higher temperature in 
the northern side is probably due to internal heating from the 
embedded YSOs. As for the southern side, the average temper- 



ature of 16 K is higher than typical temperatures for low-mass 
externally heated starless cores (~10 K, e. g., Tafalla et al. 2002, 
2004). Such higher temperatures could indicate that there are 
low-mass non-detected YSOs heating the southern side of the 
main cloud, or that there is external heating. A similar result was 
found by Li et al. (2003), who derived temperatures of ~15 K 
toward massive (with masses similar to the mass derived for the 
main cloud, see next paragraph) quiescent cores, with no signs 
of star formation, in Orion. Li et al. (2003) find that the tempera- 
tures of the massive quiescent cores can be well explained by the 
dust being heated by the external UV field from the Trapezium, 
at a distance of ~1 pc from the massive cores. In our case, the 
Bl star associated with the UCH11 region (see § 15. 3t could be 
the source of heating of the southern side of the main cloud, at 
~0.5 pc from the UCH 11 region. However, we find no temper- 
aure gradient across the south of the main cloud. Alternatively, 
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the heating could be produced by the nearby O and B stars of the 
Cygnus OB2 association (Le Duigou & Knodlseder 2002). 

5.1.2. Column density and mass 

The total mass of gas and dust derived from the dust continuum 
emission at 3.15 mm, integrated in a region including all the 
BIMA sources, is -100 Mq (assuming T d = 25 K, and /? = 1.5), 
about half of the total mass from dust at 1 .2 mm observed with 
the IRAM 30 m by Beuther et al. (2002a, 2005a), 230 Mq. We 
compared this mass with the mass derived from the high-density 
tracers. The average NH3 column density for the main cloud, 
1.45xl0 15 crrT 2 , is similar to the values obtained by Jijina et al. 
(1999) for starless regions, and the derived average N2H + col- 
umn density of 1.21xl0 13 cirT 2 agrees with the N2H + column 
densities obtained by Pirogov et al. (2003) for molecular cloud 
cores with massive stars and star clusters. The mass estimated for 
the main cloud from NH3 is ~340 Mq, assuming an abundance 
for NH3 of 10~ 8 , and the mass derived from the N2H + column 
density map is 75-190 Mq for N2H + abundances in the range 
(2-5)xl0~ 10 (Pirogov et al. 2003; Tafalla et al. 2004). Thus, al- 
though the uncertainty in the abundances is certainly affecting 
the estimation of the mass from NH3 and N2H + (indeed, there is 
strong chemical differentiation across the main cloud, see § 14.31 ), 
we obtained similar masses from the dense gas tracers and the 
dust emission. Regarding the western cloud, the derived mass 
from NH 3 and N 2 H + lies in the range 1 1^47 Mq. 
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Fig. 9. N 2 H + p-v plot for the F 1 = 0-1 hyperfine at P.A. = 1 10°, 
that is, along the western cloud and the northern side of the main 
cloud. Central position is a(72000) = 20 h 31 m 13 s .63; 5(72000) = 
+40°03'16'/7 (methanol peak), and positive position offsets are 
toward southeast. Contours start at 0. 15 Jy beam -1 , increasing in 
steps of 0.15 Jy beam -1 . 



5.1.3. Starless core candidates 

We identify three starless core candidates in the region: BIMA 3, 
BIMA 4, and the western cloud. For these three cores, the NH3 
column density is high, > 3 x 10 15 cm -2 , the rotational tem- 
perature is low, ~14 K, there is no near-infrared emission in 
the J, H, and K bands of 2MASS, no mid-infrared emission in 
the 24 /im MIPS (Multiband Imaging Photometer for SpitzeQ) 
band, and no signs of molecular outflows emanating from the 
cores (see § 15.21 for a discussion on the outflow emission in 
BIMA 4). Obviously, these cores could harbor low-mass YSOs 
(with masses < 1 Mq, as estimated from the available data) non- 
detected at millimeter and infrared wavelengths. However, the 
rotational temperature map of Fig.fTOb shows local minima asso- 
ciated with these cores, specially with BIMA 4 and the western 
cloud, supporting that the heating is external and thus that there 
is no embedded YSO yet. 

It is worth noting that, for the three cores, line widths derived 
from NH3 (and also N2FP) hyperfine fits are around 1 km s -1 
(Tabled, which are considerably larger than typical line widths 
measured toward starless cores in low-mass star-forming re- 
gions, ~0.2 km s -1 (see, e. g., Crapsi et al. 2005). For the massive 
quiescent cores studied by Li et al. (2003) the NH3 line widths 
were between 0.6 and 0.9 km s _1 , similar to what we find. Since 
the thermal component is around ~0.2 km s -1 at 15 K, the line 
widths that we found suggest that these cores are dominated by 
turbulence, as Li et al. (2003) state, or that there are important 
unresolved (5 0.05 pc) systematic motions. 

5.1 .4. Effects of chemical abundance and density 

The NH3/N2FL ratio map of Fig. [TOt l shows a clear gradient 
in the ratio across the cloud, from northwest, with low values, 
to southeast, with the highest values. In particular, there seems 

3 see http://ssc.spitzer.caltech.edu/mips/ 



to be an anticorrelation between the NH3/N2FL ratio and the 
evolutionary stage of the sources embedded in the cloud. On one 
hand, in the north of the main cloud we found most of the YSOs 
in the field, where the value of the ratio is < 100. On the other 
hand, the ratio rises up to ~300 in the southeastern side of the 
main cloud and in the western cloud, where we found starless 
core candidates. The values of the ratio in the entire region are 
consistent with those derived for low-mass star-forming regions 
by other authors (Caselli et al. 2002; Hotzel et al. 2004), who 
also find the same trend of low values (about 60-90) associated 
with YSOs, and high values (140-190) associated with starless 
cores. 

One could consider whether the high values for the 
NH3/N2FL ratio in the southeastern side of the main cloud and 
in the western cloud are just reflecting the different excitation 
conditions of NH3 and N2FL, since the critical density of NH3, 
~10 4 cm -3 (e.g., Ho 1977), is lower than the critical density 
of N2H + , ~10 5 cm -3 (e.g., Joergensen 2004), and thus in re- 
gions where the density is ~ 10 4 crrT 3 , NH3 could be thermalized 
while N2H + is faintly detected because not thermalized. We did 
a rough estimation of the density in the southeastern side of the 
main cloud from the continuum emission observed at 3.15 mm 
toward BIMA 3, adopting the mass derived in this work (see 
Table [2]) and a size of ~10", and obtained a volume density of 
~10 5 cm -3 . For this density, the transitions of both NH3 and 
N2H + are thermalized enough to be detectable. Furthermore, as 
the effects of opacity and different T sx in the cloud were already 
taken into account in the calculation of the column densities, the 
different excitation conditions of NH3 and N2H + cannot play an 
important role in the large gradient of one order of magnitude in 
the NH3/N2H + ratio across the main cloud. 

A possible explanation for the higher NH3/N2H + ratio to- 
ward starless cores may be given in terms of abundance and 
density. For long-lived (~10 6 yr) cores with moderate densi- 
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Fig. 10. a) Rotational temperature map from NH3 (1,1) and 
NH3 (2,2) obtained for the main cloud and the western cloud 
of the IRAS 20293+3952 region. Scale units are in K. b) NH 3 
column density map. Scale units are in crrT 2 . c) N2H + column 
density map convolved to an angular resolution of 7". Scale units 
are in cm" 2 , d) N(NH3)/N(N2H + ) ratio map after convolving the 
NH3 column density map from panel b to a final angular resolu- 
tion of 7". Symbols are the same as in Fig. [2] Beams are shown 
in the bottom right corner, and are 6'.'9 x 3'.'0, with P.A.= 71°5 
for panels a and b, and 7" for panels c and d. 



ties (10 4 — 10 5 cm" 3 ), there is a 'golden period' of high abun- 
dance for NH3 (Aikawa et al. 2005; S. Viti, private communi- 
cation), while N2H + keeps a constant abundance, as has been 
seen toward starless cores (e.g., Tafalla et al. 2004). This is 
consistent with the NH3/N2I-F ratio being high in the southern 
side of the main cloud, where we find starless cores with faint 
dust emission, indicating that the volume density is moderate. 
For higher densities (~10 6 cm 4 ) one may consider two mech- 
anisms to explain the low NH3/N2H + ratio. On one hand, the 



high-density cores, which have short free-fall timescales, may 
not have enough time to produce an amount of NH3 comparable 
to the southern (moderate-density) cores. On the other hand, the 
NH3 may start to freeze out onto dust grains, while N2I-F seems 
to remain unaffected up to densities of ~10 7 cm" 3 (Aikawa et al. 
2003, 2005). Note that the temperature in the main cloud is on 
average low enough so that freezing out of NH3 is possible (hot 
spots are very localized). In the main cloud the highest densities 
are reached in the northern side, where there is the strongest dust 
emission. We did a rough estimation of the average volume den- 
sity toward BIMA 1, from the derived mass and the size (at the 
5cr contour level) of BIMA 1 (Table and obtained an average 
volume density toward BIMA 1 of 5 x 10 5 cm" 3 . Consequently, 
we find low values of the ratio associated with high volume den- 
sities, and thus, with most of the YSOs in the main cloud. 



5.2. Interaction of the YSOs with the surrounding gas 

5.2.1. The UCH n region: 

In Fig. QT| we overlay the CN(l-O) emission from Beuther et 
al. (2004b) on the NH3 (1,1) emission. The morphology of the 
CN emission next to the UCH n region is following the edge of 
the dense cloud traced by NH3 (1,1), and is facing the UCH 11 re- 
gion. CN emission is often found in photon-dominated regions, 
tracing the part of the molecular cloud which is exposed to an 
intense flux of FUV photons from a nearby young OB star (see, 
e.g., Fuente et al. 1993; Simon et al. 1997; Boger & Sternberg 
2005). In addition, we found heating toward this edge of the 
cloud (Fig. [T0b), with rotational temperature progressively in- 
creasing when approaching the UCH n region. All this indicates 
that the UCH 11 region is radiatively interacting with the main 
cloud. The mechanical interaction of the compression front of 
the UCH n region with the main cloud is likely traced by the H2 
emission at 2.12 fim (Kumar et al. 2002; Beuther et al. 2004a), 
which reveals a ring-like structure around the UCH n region, of 
radius ~8", reaching the edge of the main cloud. 

5.2.2. Outflow A: 

Beuther et al. (2004a) discovered a highly collimated outflow 
in CO (2-1) (outflow A), elongated in the northeast-southwest 
direction, with mml being its driving source. The axis of the 
outflow is shown in Fig. [2] 

The velocity gradient seen in CH3OH toward mml (see § 13.21 
and Fig. |6]l is consistent with tracing material entrained by ou- 
flow A. In addition, the p-v plot shows a clump at offset posi- 
tion ~ -35", the methanol western clump, with wide redshifted 
wings (see Figs. [3] and and matching well the red lobe of 
CO (2-1) and the redshifted clump of SiO(2-l) from Beuther 
et al. (2004b). The methanol western clump appears in the axis 
of outflow A, just downstream of a weak dust filament that con- 
nects IRS 1 and IRS 2 (see Figs. Q] and [2^). This suggests that the 
interaction of outflow A with the filament sweeps-up and chemi- 
cally alters the material of the filament, originating the methanol 
western clump. 

The red lobe of outflow A, which starts at mml, propagates 
through the cavity apparent in the N2H + and NH3 zero-order 
moment maps (see Fig. The rotational temperature and col- 
umn density maps of Fig. [10] show that the walls of this cavity 
are characterized by remarkably high temperatures and low col- 
umn densities. In addition, the NH3 gas of the walls of the cavity 
is redshifted (as the lobe of outflow A), and shows line broad- 
ening. Thus, this cavity could have been excavated by outflow 
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Fig. 11. Zero-order moment map of the CN (1-0) emission (con- 
tours) from Beuther et al. (2004b) in the IRAS 20293+3952 
region, overlaid on the NH3(1,1) emission shown in Fig. |2j; 
(grey scale). CN(l-O) emission has been integrated from 4 to 
12 km s , and contours start at 15% of the peak intensity, 
273 Jy beam -1 km s , increasing in steps of 15%. Symbols are 
the same as in Fig. [2] Beams are 6'.'9 x 3'.'0, with P.A.=72° for 
NH 3 (1,1), and 1'.'46 x 1'.'21, with PA.=32° for CN (1-0). 



A, and the high temperatures could arise from the shock inter- 
action of outflow A with the walls of the cavity. This has been 
observed toward other outflows driven by high-mass stars (e. g., 
AFGL 5142: Zhang et al. 2002). 

5.2.3. Outflow B: 

Outflow B was first identified by Beuther et al. (2004a) in 
SiO (2-1), which shows a very strong blue lobe elongated in the 
northwest-southeast direction, and in CO (2-1), with faint high- 
velocity blueshifted emission. CH3OH emission from outflow B 
is morphologically very similar to the SiO (2-1) blue lobe and 
also splits up into two lobes in a fork-like structure at the po- 
sition of BIMA 4. A scenario in which outflow B is interacting 
with a starless core, BIMA 4, naturally explains the partial de- 
flection of the outflow seen in the SiO and CH3OH emission. 
There are several pieces of evidence that support this scenario. 
First, the CH3OH emission splits up and has a change in kine- 
matics when it reaches BIMA 4 (§ 13.2b . Second, at this posi- 
tion we found some heating and the NH3/N2I-F abundance ratio 
shows lower values than those in the southern part of BIMA 4 
(Fig. [TOk.d; temperature enhancement and chemical variations 
are signatures of shocked molecular gas, e. g., LI 157: Bachiller 
& Perez Gutierrez 1997). Third, the N2I-F line broadening along 
outflow B (which is also observed in NH3) indicates that the out- 
flow is already interacting with the molecular cloud (Fig. [8}?). 
Finally, downstream of BIMA 4, there is heating (see Fig.fTUb) 
associated with the H2 knot c of Kumar et al. (2002). The partial 
deflection of an outflow due to the interaction with a dense qui- 
escent clump has been found toward other star-forming regions 
(e. g., IRAS 21391+5802, Beltran et al. 2002). Interaction of an 
outflow with a dense flattened NH3 core has been reported too 
toward NGC 2024 by Ho et al. (1993), where the outflow sweeps 
the material off the surface of the NH3 core. Similarly, outflow 
B interacting with the BIMA 4 core could be also sweeping the 



material off the surface of the core, producing all the observa- 
tional features described above. 



In addition, the properties of the methanol eastern clump (lo- 
cated 25" to the northeast of BIMA 4), which is found at ambi- 
ent velocities, and with narrow lines (around 1 km s , Fig. [3}, 
could be a consequence of the illumination by the UV radia- 
tion coming from the interaction of outflow B with BIMA 4. 
The illumination of the UV radiation from shocks has been pro- 
posed as the mechanism of enhancement of the emission of some 
species in clumps ahead of shocks, in particular of CH3OH and 
NH 3 (see e.g., Torrelles et al. 1992b; Girart et al. 1994, 2002). 
With the present observations, we were not very sensitive to the 
NH3 emission toward the methanol eastern clump because the 
clump is located beyond the VLA primary beam FWHM (see 
Fig.[2j;,d). Typically, such illuminated clumps in low-mass star- 
forming regions have sizes of 0.05-0.1 pc, and are located at 
distances to the shock of the order of 0.1 pc (Girart et al. 1998). 
This is consistent with the methanol eastern clump, which has a 
size of ~0. 15 pc, and is located about 0.2 pc from BIMA 4. 

It is interesting to note that the morphological features of the 
dense gas around outflow B are different from outflow A. While 
outflow A has already cleared up the dense molecular gas, cre- 
ating a well defined cavity, outflow B is probably in the embri- 
onary phase of creating the cavity. This suggests that possibly 
outflow B is younger or less energetic than outflow A, in agree- 
ment with the measurements of Beuhter et al. (2004a). It is worth 
noting that the line width of the dense gas associated with out- 
flow B is significantly higher than the thermal line width for both 
NH3 and N2H + (~ 1.8 km s compared to 0.2 km s~'). A similar 
result is found by Wang et al. (2006) toward massive protostel- 
lar cores mapped in NH3 with high angular resolution. Given 
that we do not observe systematic motions in the dense gas as- 
sociated with outflow B, we suggest that the large line width in 
the dense gas is due to turbulence injected by the passage of the 
outflow (see Fig. [8})). 

Regarding the driving source of outflow B, a possible candi- 
date is mml, as Beuther et al. (2004b) proposed. If this was the 
case, mml would be a binary system of jet sources, since outflow 
A is clearly associated with mml. Another possibility would be 
that BIMA 4 is the driving source of outflow B. However, the 
kinematics of the gas as traced by CH3OH show that there is no 
clear symmetry with respect to BIMA 4 (see Fig. [5^: BIMA 4 
is at an offset position ~20"). Another candidate of being the 
driving source of outflow B is the Northern Warm Spot, which is 
aligned with outflow B, and is located at the begining of the blue 
lobe of CH3OH. Furthermore, we found a clear symmetry with 
respect to the Northern Warm Spot in the p-v plot of Fig.|5h (see 
§ 13.21 ), strongly suggesting that the Northern Warm Spot could 
be the driving source of outflow B. In § |5.3| we estimate an asso- 
ciated mass of ~0.7 Mq for the Northern Warm Spot. 

Finally, the interaction of outflow B with the BIMA 4 core 
leads us to speculate that outflow B could be triggering the col- 
lapse in this core, as has been proposed in other regions (e. g., 
Yokogawa et al. 2003). This would draw a scenario in which a 
YSO in the north of the main cloud is 'responsible' for star for- 
mation in the south, but should be further investigated by study- 
ing in more detail the morphology and kinematics in the core. 
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Fig. 12. (J - H), (H - K) diagram for the 2MASS sources in the 
IRAS 20293+3952 region lying within the BIMA primary beam. 
We selected those sources, IRS 2 to IRS 6, with (H-K) > 2, 
(J—H) < 3, and spatially coincident with the dense gas emission. 



very low column densities, as expected for gas after the passage 
of an ionization front (see, e. g., Dyson & Williams 1997). 



IRS 2: About 10" to the northwest of the UCHu region, there 
is the second brightest infrared source in the field, IRS 2, which 
shows much more infrared excess than IRS 1, as already stated 
by Kumar et al. (2002). We could not determine a temperature 
and column density from NH3 toward IRS 2 due to the low 
S/N of the NH3 emission, but we detected continuum emission 
at 3.15 mm, and we estimated an associated mass of 0.6 Mq. 
The high near-infrared flux of IRS 2 suggests that this is an 
intermediate-mass YSO. 



IRS 3: IRS 3 is the only infrared source with infrared excess 
falling inside the 2cr contour level of BIMA 1 . The high angular 
resolution observations of Beuther et al. (2004b) reveal a weak 
extension of mml toward the southwest reaching the position 
of IRS 3. From our data, the dense gas emission was too weak 
to derive a temperature and column density toward this source. 
Thus, IRS 3 shows strong infrared emission and weak millimeter 
emission, indicating that IRS 3 is likely a YSO which is clearing 
up the material in its surroundings. 



5.3. YSOs in the region 

5.3.1 . Selection of 2MASS sources associated with the 
region 

In order to study the different sources associated with the dense 
gas found around the UCH n region, we extracted a sample of 
stars within the BIMA primary beam from the 2MASS Point 
Source Catalogue (PSC) (Skrutskie et al. 2006), and plotted a 
(/ - H), (H - K) diagram (Fig.[T2l. The total amount of infrared 
sources inside the BIMA primary beam is 43. The color-color 
diagram shows that there is a bulk of infrared stars at low-to- 
moderate values of the color indices, which occupy the posi- 
tion of stars with no infrared excess, and a group of stars with 
high (H - K) and low (J - H) colors. We took the criteria of 
(H-K) > 2, (J - H) < 3 and spatial coincidence with emis- 
sion detected in this work to select those 2MASS sources that 
are possibly (but not necessarily) associated with the dense gas 
around the UCH 11 region. The selected 2MASS sources are la- 
beled as IRS 1, IRS 2 (following the nomenclature of Kumar 
et al. 2002), and IRS 3 to IRS 6 (increasing RA). These in- 
frared sources (with the exception of those that seem to be back- 
ground/foreground objects, see below) are listed in Table [4] to- 
gether with the compact millimeter sources from Beuther et al. 
(2004b), and the sources found in this work. For each source we 
show the main properties derived in this work. 

5.3.2. Individual YSOs 

IRS 1 and the UCHu region: From the 2. 12 /mi image provided 
by Kumar et al. (2002), we find that IRS 1 is probably a binary 
system separated ~3", the fainter component being associated 
with the UCH 11 region. The UCH 11 region has a flux density of 
7.6 mJy beam -1 at 3.6 cm (Sridharan et al. 2002; Beuther et al. 
2004a), and a deconvolved size of 4'.'5 (0.04 pc). This corrre- 
sponds to an ionizing flux of ~1.5xl0 45 s _1 , typical of stars of 
spectral type Bl (Panagia 1973). Toward the UCHu region, we 
found the highest rotational temperature in the main cloud, and 



IRS 5: IRS 5 falls inside the 7cr contour level of the dust ridge 
(see Fig. Q]), in its eastern edge. We found heating associated 
with IRS 5 (§ 14.21 and Table |3j, indicating that this source is 
physically associated with the main cloud. In addition, the NH3 
column density is low (Fig. [T0b), and thus IRS 5 could be a 
protostar in the process of clearing up the surrounding material. 
Alternatively, the low NH3 column density could be due to de- 
pletion of NH3 onto dust grains. 



mm1 : mml is likely the main contribution to the emission of 
BIMA 1, since it is the strongest compact millimeter source 
in the region, and is the driving source of outflow A, which is 
highly collimated (Beuther et al. 2004a) and elongated in the 
same direction as BIMA 1 (Fig. [TJ. In addition, the infrared 
emission at 2.12 /mi associated with mml is faint and arises 
mainly from shocked gas (see Fig. 2 of Kumar et al. 2002). The 
rotational temperature toward mml is 24 + 3 K (Table [?}, and 
the NH3 column density is rather low, possibly due to depletion 
of NH3 onto dust grains. From the flux density of Beuther et al. 
(2004b) at 2.6 mm, we obtain a mass for mml of 4.0 Mq, as- 
suming f3 = 1 (value estimated by Beuther et al. 2004b), and 
Ta — 32 K (estimated from the rotational temperature derived in 
this work). All this indicates that mml, in contrast with IRS 3, is 
deeply embedded in a massive and not very hot envelope, sug- 
gesting that it is an intermediate/high-mass protostar. 



mm2 and mm3: Beuther et al. (2004a) detect two compact mil- 
limeter sources inside the 5cr contour level of BIMA 2, mm2 
and mm3, and they suggest that each one is driving a molecular 
high-velocity CO (2-1) outflow (outflows C and D), which in- 
dicates that mm2 and mm3 are truly protostars, and not heated 
clumps of dust without any star yet. Taking the flux density from 
Beuther et al. (2004b) at 2.6 mm, a dust emissivity index j3 — 1, 
and 7d = 30 K and 7^ = 23 K for mm2 and mm3 respectively 
(7d derived from this work, as in the case of mml), we estimate 
associated masses of 1.3 Mq for mm2 and 0.8 Mq for mm3. 
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Table 4. List of sources in the 120293 region 



Position" Av(NH 3 ) c iV(NH 3 ) T m if c M c 





a(J2000) 


5(J2000) 


Outf b 


(km s- 1 ) 


(10 15 cm" 2 ) 


Ratio d 


(K) 


(mjy beam ') 


(Mq) 


YSOs 




















IRS 1 


20:31:11.26 


+40:03:07.6 


NN 


_ 


_ 


_ 


_ 


3.9 


0.6 


UCHu 


20:31:11.22 


+40:03:11.0 


NN 


2.4 ±0.3 


0.14 


57 


38 ± 15 


4.0 


0.4 


IRS 2 


20:31:10.32 


+40:03:16.5 


NN 


_ 




_ 


_ 


3.5 


0.5 


IRS 3 


20:31:12.48 


+40:03:20.0 


YN 


_ 


_ 


_ 


_ 


9.5 


1.4 


IRS 5 


20:31:13.41 


+40:03:13.7 


?N 


1.53 ±0.06 


0.59 


42 


24 ±3 


5.2 


1.2 


mml 


20:31:12.9 


+40:03:22.9 


YY 


1.37 ±0.12 


1.1 


37 


24 ±3 


28.0 s 


4.0 


mm2 


20:31:12.01 


+40:03:12.2 


YN 


1.10 ±0.07 


0.88 


46 


23 ±4 


8.68 


1.3 


mm3 


20:31:12.41 


+40:03:10.5 


YN 


1.23 ±0.13 


0.66 


57 


19 ± 1 


3.6 s 


0.8 


NWS h 


20:31:13.23 


+40:03:19.5 


NY 


2.00 ±0.11 


1.8 


56 


29 ±6 


3.7 


0.7 


SWS h 


20:31:13.26 


+40:02:58.9 


YN 


<1.03±0.04 


1.3 


150 


20 ±2 


<0.28 


< 0.9 


Starless Cores 




















BIMA3 


20:31:13.92 


+40:03:00.9 


NN 


1.18 ±0.06 


5.4 


290 


14 ±2 


4.0 


2.2 


BIMA4 


20:31:14.40 


+40:03:07.4 


YY 


1.34 ±0.06 


3.8 


180 


14 ±2 


3.5 


2.0 


WCloud h 


20:31:10.71 


+40:03:30.3 


NN 


<1.03±0.02 


2.8 


280 


15 ±3 


<0.28 


< 1.3 



a Positions correspond to those in the 2MASS PSC for IRS 1 to IRS 5; the centimeter source given by Beuther et al. (2002b) for the UCHu 
region; the positions given by Beuther et al. (2004b) for mml to mm3; the peak in the rotational temperature map of FigllOfa for the NWS and 
the SWS, and the position of the peak in the NH3 column density map of FigllOb for the starless cores. 

b High-velocity outflow associated (Y=Yes; N=No): the first character refers to CO from Beuther et al. (2004a), the second character refers to 
high-velocity CH3OH from this work. 

c Line width derived from the fits to the hyperfine structure of NH 3 (1,1) (see § !4.2t . 

d Ratio of JV(NH 3 ) over iV(N 2 H + ). 

c Intensities at 3. 15 mm from this work. The rms noise is 0.7 mly beairT 1 , and the upper limits for non-detections have been set to Act. 

f Associated mass (essentially circumstellar) from the mm continuum emission listed in column (9), assuming a dust emissivity index /3 = 1, 

and estimating the kinetic temperature from the rotational temperature listed in column (8), following the expression of Tafalla et al. (2004, 

see § 1531 

8 For mml, mm2, and mm3 the intensity corresponds to the flux density (in mjy) given by Beuther et al. (2004b). 
h NWS: Northern Warm Spot; SWS: Southern Warm Spot; WCloud: western cloud (see Fig.|2p). 



The Northern Warm Spot: The Northern Warm Spot is located 
about 5" (10000 AU) to the southeast of mml, and is near a local 
maximum in NH3 and N2H + column densities (see Fig.fTOb and 
c). From the millimeter continuum emission we estimated an as- 
sociated mass of ~0.7 Mq (assuming /3 — 1 and Ta - 44 K). 
Very close to the warm spot, l'/3 to the east, there is faint in- 
frared continuum emission at 2.12 /mi, which is not detected in 
the J and H band images of 2MASS, indicating that the source 
has strong infrared excess. Since there is no H2 line emission at 
2.12 yum tracing shocked gas (see Fig. 2 of Kumar et al. 2002), 
the infrared emission is most likely tracing an embedded proto- 
star, and not a condensation heated up by the impact of outflow 
B. These properties support that the Northern Warm Spot is the 
driving source of outflow B. 



The Southern Warm Spot: Similar to the NWS, the Southern 
Warm Spot is found ~10" to the southwest of BIMA 3 (see 
jj4.21 >. At this position, the temperature enhancement is very sig- 
nificant (see Table |3), and we also found faint infrared emis- 
sion associated, present only in the K filter of 2MASS, and with 
no H2 line emission associated. The emission from the con- 
tinuum at 3.15 mm is below 4cr, implying an associated mass 
< 0.9 Mq (assuming = 25 K, estimated as in the case of 
mml). Furthermore, there is high-velocity blueshifted CO (2- 
1) emission spatially coincident with the Southern Warm Spot 
(Beuther et al. 2004a). All this indicates that possibly the heat- 
ing in the Southern Warm Spot is due to an embedded low-mass 
YSO. 



Background/foreground objects: In addition to IRS 5, we find 
two 2MASS sources with infrared excess, IRS 4 and IRS 6, spa- 
tially coincident with the dust ridge. IRS 6 is the source with the 
highest infrared excess among the 2MASS sources in the main 
cloud, suggesting that it is embedded in gas and dust. The faint 
source IRS 4 is located only 4" to the southwest of IRS 6, and 
also shows infrared excess. However, the association of IRS 4 
and IRS 6 with the dense gas is not clear from the rotational tem- 
perature and NH3 column density maps (Fig. [TOk.b), and thus 
these two infrared sources could be highly extincted background 
stars or foreground YSOs belonging to the cluster. 



Summary: As seen above, there is a variety of YSOs in the re- 
gion. The most massive sources seem to be the UCH n region, 
IRS 1, and IRS 2. To the east of the massive sources, there 
is the main cloud of dense gas, where we find seven sources 
(IRS 3, IRS 5, mml, mm2, mm3, the Northern Warm Spot 
and the Southern Warm Spot), which seem to be (proto)stars 
coming from the same natal cloud. In addition, we found three 
cores with characteristics similar to starless cores. The sources 
seem to be in different evolutionary stages, with the infrared 
sources in the more advanced phases, the millimeter sources 
with no infrared emission likely being in an earlier stage, and 
the youngest sources having the properties of starless cores. We 
do not find different evolutionary stages only between the low- 
mass sources and the high-mass sources (which are expected be- 
cause the high-mass sources evolve more rapidly to the main 
sequence than the low-mass sources), but among the low-mass 
sources as well, for which we expect similar rates of evolu- 
tion to the main sequence. For example, IRS 3 and IRS 5 are 
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YSOs strongly emitting in the infrared but with faint millime- 
ter emission, and do not seem to be driving outflows, sugges- 
tive of being at the end of the accretion phase. Contrary to this, 
mm2 and mm3 are YSOs with the millimeter emission much 
stronger than the infrared emission, and they are the driving 
sources of collimated outflows, indicative of being in the main 
accretion phase. The masses of IRS 3 and IRS 5 are difficult to 
estimate, but comparing with the infrared emission of the known 
intermediate/high-mass sources one may classify them as low- 
mass YSOs. The masses of mm2 and mm3, estimated from the 
millimeter emission, are around 1 Mq. Thus, these low-mass 
YSOs (IRS 3/IRS 5, and mm2/mm3) seem to be in different 
evolutionary stages, and we conclude that stars are not forming 
simultaneously in this cluster environment. Rather, there may 
be different generations, as found toward other star-forming re- 
gions (e. g., NGC 6334: Beuther et al. 2005b; S235A-B: Felli et 
al. 2006; L1551: Moriarty-Schieven et al. 2006). 



5.4. Spatial distribution of the YSOs in the region 

In this section we consider whether interaction between the dif- 
ferent sources is important in the determination of the spatial 
distribution of the YSOs. We find that star formation is localized 
in the north of the main cloud, where there are around six YSOs, 
while in the south we find sources with properties of starless 
cores, and only one YSO. 

We consider first whether the UCH n is responsible for such a 
spatial distribution. In section l3T2l we discussed some evidence of 
interaction of the UCH n region with the edges of the main cloud, 
mainly with BIMA 2. This suggests that mm2 and mm3 could 
have been triggered by the UCH n region. If this was the case, 
one would expect to find some evidence of a compression front 
expanding away from the UCH n region. In the p-v plot of Fig. [9] 
the N2H + emission shows the shape of an incomplete ring, which 
could be interpreted as an incomplete expanding shell. However, 
the velocity field as seen from N2H + in Fig.[8]3 does not show any 
radial symmetry with respect to the UCH n region, and the pres- 
ence of outflow B makes difficult to disentangle the advance of 
any compression front. In addition, the formation of the YSOs in 
the northeast of the main cloud is not likely due to triggering by 
the UCH ii region, since they are located farther away than mm2 
and mm3 and are not in earlier evolutionary stages. Therefore, 
triggering by the UCH n region does not seem to be the dominant 
agent causing star formation in this clustered environment. 

Another possibility for the formation of the YSOs in the 
north of the main cloud could be a merging of the main cloud 
with the western cloud. In Fig. [9] we found that the western 
cloud has an extension intersecting the main cloud, seen at ve- 
locities > 7 km s _1 in the channel maps (Pig. [7j. Merging of two 
clouds has been proposed as a mechanism to trigger star forma- 
tion in other regions, (Wiseman & Ho 1996; Girart et al. 1997; 
Looney et al. 2006; Peretto et al. 2006). Note however that no 
further evidence of such a merging (heating, line broadening, 
two clear velocity components across the main cloud) is seen as- 
sociated with the extension of the western cloud, and thus the 
merging scenario is not fully consistent with our data. 

We therefore conclude that the spatial distribution of the 
YSOs in this cluster environment could simply reflect the initial 
conditions of the cloud: if the main cloud was originally much 
denser in the north than in the south, we would possibly observe 
a similar situation to what we have found. 



6. Conclusions 

We observed with the BIMA and VLA arrays the continuum 
emission at 3.15 mm, and the N2H + , NH3, and CH3OH emis- 
sion toward IRAS 20293+3952, a region in which star forma- 
tion is taking place in a closely-packed environment. Our main 
conclusions can be summarized as follows: 

1. The dense gas traced by N2H + and NH3 shows two differ- 
ent clouds, one to the east of the UCH 11 region (main cloud), 
of ~0.5 pc of size and ~250 Mq, and another cloud to the 
northwest (western cloud), of ~0.15 pc and ~30 Mq, and 
redshifted with respect to the main cloud. The dust emission 
reveals two strong components in the northern side of the 
main cloud, BIMA 1 and BIMA 2, and two fainter compo- 
nents in the southern side, BIMA 3 and BIMA 4, together 
with extended dust emission forming a common envelope. 
Regarding the CH3OH, we found strong emission in a fork- 
like structure associated with outflow B from Beuther et al. 
(2004a), as well as two CH3OH clumps associated with out- 
flow A. 

2. We found that the rotational temperature is higher in the 
northern side of the main cloud, around 22 K, than in the 
southern side, around 16 K. In contrast, the NH3 column den- 
sity has the highest values in the south of the main cloud. The 
N2H + column density distribution resembles the dust emis- 
sion, strong in the northern side of the main cloud. We found 
three local temperature enhancements which seem to be as- 
sociated with embedded YSOs, one of them associated with 
a 2MASS source, and the other two, the Northern Warm Spot 
and the Southern Warm Spot, associated with faint contin- 
uum infrared emission at 2.12 fim. 

3. There is strong chemical differentiation in the region. In par- 
ticular, we found low values of the NH3/N2H + ratio, ~50, 
associated with YSOs, and high values, up to 300, asso- 
ciated with starless cores. This is consistent with NH3 be- 
ing enhanced with respect to N2H + at moderate densities 
(< 10 5 cm -3 ), while at densities around ~10 6 cm 4 the NH3 
enhancement may be less important (due to shorther life- 
times of the cores), and additionally the effects of NH3 de- 
pletion may start to play a role. 

4. We identified three cores in the NH3 column density map, 
with low temperatures and with no infrared emission asso- 
ciated, similar to starless cores in low-mass star-forming re- 
gions, but with temperatures slightly higher, ~15 K, possibly 
due to external heating from the nearby OB stars. These are 
associated with BIMA 3, BIMA 4, and the western cloud. 

5. Interaction between the different sources in the region is im- 
portant. First, the UCH 11 region is interacting with the main 
cloud heating and enhancing the CN(l-O) emission in the 
edge of the main cloud that is facing the UCH 11 region. 
Second, one of the outflows in the region, outflow A, seems 
to be excavating a cavity and heating its walls. Third, an- 
other outflow in the region, outflow B, is interacting with the 
BIMA 4 starless core, likely producing the deflection of the 
outflow seen at this position. Such a shock could be illumi- 
nating a clump with narrow CH3OH lines located ~0.2 pc 
to the northeast of the shock. We propose that the Northern 
Warm Spot may contain the driving source of outflow B. The 
properties of the molecular gas around the outflows suggest 
that outflow A is more energetic and/or older than outflow B. 

6. There are about eight YSOs in the dense gas near the UCH 11 
region. For the YSOs for which we could estimate the associ- 
ated circumstellar mass, we found values ranging from < 0.4 
to ~4 Mq. In addition, the YSOs seem to be in different 
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evolutionary stages, even if we consider only the low-mass 
sources. Thus, stars do not seem to form simultaneously in 
this cluster environment. 
7. While we cannot discard that interaction between the differ- 
ent sources may have triggered star formation in particular 
cases, triggering cannot explain the overall spatial distribu- 
tion of the YSOs, suggesting that this may be essentially de- 
termined by the initial conditions in the cloud. 
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